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ABSTRACT 

Unitarity for the first row of the quark mixing matrix appears to be fully satisfied 
by the value of V^s obtained from recent new measurements of kaon semileptonic 
decay rates. 



1 Introduction 



What is today called quark mixing was first introduced in 1963 by Cabibbo jl] as an 
angle characterizing the relative strength of strangeness changing and strangeness 
conserving weak amplitudes. In modern language we write the quark weak current 
as Jq, = U7a(l — 75)VD, where U, D are the charge 2/3, —1/3 quark vectors and 
V is a unitary matrix |2], often called the CKM matrix. Whether V is unitary can 
be checked experimentally. 

For no good reason at all it became fashionable, in the mid-late nineties, to 
expect the CKM-matrix not to be unitary, providing great impetus to the study of 
B mesons. The beautiful experiments at KEK and SLAG have dashed such expecta- 
tion, verifying unitarity to (9(10%). By far the best check comes however from kaon 
physics. The so-called first row unitarity requires |Kidp + + |Ki6p = 1 which, 
because iKitp ~ 10~^ jS], is equivalent to iV^dp + \Vus\^ = 1- Semileptonic kaon 
decays offer possibly the cleanest way to obtain an accurate value of the Cabibbo 
angle or better \Vus\- Since i^' — ^ vr is a O^^O^ transition, only the vector part of 
the weak current has a non vanishing contribution. Such processes are protected by 
the AdemoUo-Gatto theorem against SU(3) breaking corrections to lowest order in 
ms—rria. Recent advances in lattice gauge techniques have allowed the computing of 
the pseudoscalar decay constants /^r and As a consequence, the K ^ jiu partial 
decay rate provides a new handle for the determination of \Vus\- 

2 Experimental situation prior to 2004 

For several years the Particle Data Group (PDG) has reported values of \Vud\ and 
\Vus\ in slight disagreement with unitarity. The 2004 version [7j, gives =0.2200 
±0.0026 and |Kd|=0.9738±0.0005 from which the sum of the squares is 0.9966 
±0.0015 which verifies unitarity to (9(0.33%) albeit with a deviation of ~2cr. This 
result is however based on the so called PDG fit. Examining inputs used for the fit 
and the output of the fit, one notices a few peculiarities: 

1. Direct measurements of the partial rates and/or branching ratios for the rele- 
vant decays, K neu, TTfiu, are quite old, with accuracies of (9(5%) 

2. By using ratios of partial rates, even far removed from the quantities of interest, 
the fit comes out with errors as low as 0.7%, this being the case for i^'^— >-7r^e^i/ 

3. Radiative contribution to the input values are essentially unknown. This in- 
validates the constraints used in the fit as well as the usefulness of the resulting 
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BR values for obtaining IV^^I 

4. There are large differences in the form factor slopes for neutral and charged 
kaons. The ratio i^esZ-ft'/iS which is a unique function of the form factor slope 
parameters is quite inconsistent between charge and neutral kaons. The most 
recent results, from KEK jl] and Istra+ OIHl, clearly indicate a smaller ratio. 

5. The charged kaon lifetime, given by PDG with an accuracy of 0.2% is in fact 
an average or fit of values covering an interval of 1.5%. While PDG might 
wish to be democratic, some of those values must be wrong. 

Both central values and errors given by PDG for branching ratios are therefore 
suspect and anyway not usable, because of ignorance about inclusion of radiation, 
to obtain \Vus\ to better than several per cent. 

3 How to get to \Vus\ 

In principle the connection between partial decay rates and \Vus\ is as simple as: 

T{K 7r£iZ/) = TV IK/ y" W\'^dEidE2 x C (1) 

where contains well known things like the Fermi constants, some tt's and two's, 
Clebsches etc. There is a phase space integral and some corrections C. We begin 
with things that follow from measurement. Partial width are always obtained from 
F = BR/r, i.e. we must measure both branching fractions and lifetimes. The phase 
space integral is a function of masses, quite well known and of the shape of the form 
factor. The /+ form factor is defined by: {7i\J^\K) = f+{t = 0) x f+{t){P + p)^ 
where = 1 + X'jtlm+ (A'_|_/2)(t/m)^ + .. or some other parameterization, with 

/+(0) = 1. P, p are the initial and final four momenta, m is the pion mass and is 
the vector part of the quark current. The form factor parameters can and must be 
measured for the vector form factor /+ and also for the scalar form factor /q. In the 
following we retain the three parameters: A'_^,A'| defined above and the analogous 
parameter Aq. 

Amongst the corrections, radiative effects are quite important. Since we 
cannot turn off electromagnetism, we must include in the measurement some well 
defined fraction of the radiative decays. Ideally we want to measure the rate or BR 
for decays to vrez/ etc., plus any number of 7's. This rate is finite and calculable. A 
very important point, being only recently fully recognized, is that radiative effects 
must be included in the acceptance calculations, in practice in the decay generators 
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used for the Monte Carlo simulations. Otherwise counting mistakes of several per 
cent are inevitable. 

There is more to the corrections. There are radiative corrections, both 
short distance and low energy corrections to the phase space density both with and 
without external radiation. There are isospin breaking correction to account for 
7T^-7i^ and K^-K^ (or u-d) differences. Finally pions and kaons are not quite the 
same objects and SU(?>)f breaking (or s-d difference induced) corrections must be 
included. This last corrections is the factor /+(0), mentioned above in the definition 
of the form factor. We can note that there is a unique value f^°{Q) for both /+ 
and /o in neutral kaon decays and another for charged kaons. One ends up with 
something like:^ 

r = ^Ew l/f (o)p cl [1 + + 51] (2) 

/j(A', A", Ao) is the integral of the phase space density, factoring out and without 
radiative corrections. Short distance radiative corrections are in the Sirlin term S'ew 
jS]. In addition long distance radiative corrections OCD] for form factor, phase space 
density and SU(2) breaking are included as and A/j(A). 

4 New indications 

The first evidence that not all was correct with the PDG information came from 
E865 at BNL, who pubhshed in 2003 a fully inclusive value for K'^^ decays [TT| of 
BR{K+ 7r0e+i/(7)) = (5.13±0.10)%, to be compared to (4.87±0.06)%. The differ- 
ence is not very large and moreover E865 actually measures the ratio e3 / (7r2 + /i3 + 
7r3) and relies on PDG for knowledge of the denominator. Early in 2004 KLOE ^2] 
reported on V{Ks TT^e^v {'-))), finding r=(7.92±0.12) x 10^ s"\ to be compared to 
(7.50 ±0.08) X 10^ s"^ for Kl- This resuh does not rely in any way on external infor- 
mation but uses the new, precise KLOE result for r(7r+7r^) /r(7r'^7r°) The re- 
sults implies |K.|=0.2272±0.0029 while from KLOE |Ks|=0.2254±0.0026, both val- 
ues in very good agreement with unitarity which requires a value of 0.2265±0.0022. 

5 Doing it correctly 

Past experience proves the dangers of using too many intermediate steps and infor- 
mation from different experiments. The best way to obtain partial rate measure- 
ments is to start with a tagged kaon sample and count decays to the appropriate 

^The numerical factor 768=2^x3, also appears sometimes as 128, other times as 192. They are 
all equally valid and the difference is made up by the phase space integrals Li. 
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channels. Accurate lifetimes are still lacking, especially for i^^j^-mesons but also pos- 
sibly for charged kaons. Finally, to extract \Vus\ the form factor parameters must 
be measured, values calculated in chiral perturbation not being reliable. 

An alternate way consists of measuring all major ratios Fj/F, where Fj is 
the partial width for channel i and F is the partial width for a major channel, not 
included in the previous set. Then the constraint F + X^Tj = 1 — e allows solving for 
all branching ratios Fj/F and F/F. For the Kl case, choosing the 4 channels Treu, 
TTfiu, 7r+7r~7r° and vr'^Tr^Tr'^ we have 0.0035 < e< 0.0036 i.e. one can reach in principle 
an accuracy of ~0.0001 for Kes and Kf^^ or ~0.03%. 

The first method is ideally employed at a ^-factory where the production 
of the KsKl pair in 0-decay provides a tagged, monochromatic i^^-beam of known 
flux. This is the choice of KLOE jHj which has obtained preliminary results dis- 
cussed below. An advantage of working with a tagged monochromatic Kl beam is 
that of performing also lifetime measurements. The latter can be done by observing 
the time dependence of decays to a particular channel. Even better, the lifetime can 
be obtained by just counting all decays in a time, or space, interval. 

Final results on decays have been given by the KTeV collaboration |15j . 
In their very complete study of semileptonic Kl decays, KTeV gives very accurate 
branching ratios and new precise values for the form factor parameters. The Kf,^ 
branching fraction is found to be 0.4067±0.0011, strikingly different from the PDG 
fit value of 0.3878±0.0027. They have also clearly established the necessity of a 
curvature term in the vector form factor, quantitatively A'l = (3.20 ± 0.69) x 10^^. 
The linear term is however smaller,^ A'_,_ = (20.64 ± 1.75) x 10"'^. This results in a 
1% reduction of the phase space integral corresponding to an increase of 0.5% for 
\Vus\- This reinforces the warning that comparison of results obtained with different 
parameter values are meaningless. 

6 The KLOE results 

Another advantage of the monochromatic tagged Kl beam is the unambiguous 
closure of the kinematics. This however is not quite so important for just counting 
the number of different decay modes. There are two main cases to distinguish: 
1) all neutral particles: Tr^Tr^Tr*^ is ~99.8% of the total; 2) two charged particles 
final state, dominated by TT^e^u, n^fi^u and 7r"'"7r~7r°. The three contribution are 
quite well distinguished by the variable A^^ = -Emiss — |Pmiss|, computed using the 
muon and the pion masses for the charged particles and choosing the lowest value. 
^Partial indications were given in 0. 
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Fig. 1 illustrates the channel separation as well as the agreement with Monte Carlo 
simulation which fully accounts for radiation in all channels. 
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Figure 1: Distribution of the variable A^^ = E^igg — |pmiss|; s^e text, for decays 
to two charged particles. Only a few per cent of the data is plotted. 

The validity of this procedure is confirmed by using particle identification by time 
of flight and range-di?/da; in the calorimeter. Fig. 2 shows an example of the Ag^r 
distribution for events with an identified electron compared to MC prediction. 

9 
I 




Figure 2: Distribution of the variable Ag^, see text, for Kl decays with identified 
electrons. 



The final results for the branching ratios are: 

BR{Kl TT^e^i/) = 0.3985 ± 0.0035 
BR{Kl TT^i^^iy) = 0.2702 ± 0.0025 



(3) 

(4) 
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BR{Kl 7r°7r°7r°) = 0.2010 ± 0.0022 (5) 

BR{KL^n+n-n^) = 0.1268 ±0.0011 (6) 

The above results are obtained from the study of ~13 milhon Kl decays. From the 
number of observed decays in the fiducial volume we obtain: 

t{Kl) = 51.35 ± 0.05,tat ± 0.26syst ns (7) 

An additional 40 million decays are used to obtain the various efficiencies. Details 
of the determination of BIi{KL—>-TT^7[^TT^) are beyond the present discussion. 

The entire data sample is used to measure the lifetime. Fig. 3 shows the 
time dependence of Kl—^tt^tt^tt^ decays in the KLOE drift chamber and the region 
used in the lifetime fit. Regeneration in the beam pipe and inner chamber wall is 
visible. 
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Figure 3: Distribution of the decay proper time for KL—>-TT^7r^7[^. 

We find: 

TK^ = 51.15 ±0.20,tat± 0.40 syst ns [O) 

These preliminary values are quite consistent with previous measurements [7j and 
together give a new improved lifetime: t{Kl) = 51.35 ± 0.2 ns. 



7 IK.I 

Before getting to \Vus\, obtained from inverting eq. El we wish to remark that ex- 
perimental errors are getting smaller than 1%. Specifically the lifetime contributes 
an error of 0.2%, the BR values about 0.15% and the form factor parameters some 
0.3%. Knowledge about the form factor normalization f^° remains relatively poor. 
Leutwyler and Roos gave in 1984 [in] the estimate f^ =0.961±0.008 which corre- 
sponds to a 0.83% uncertainty. This result is confirmed by by lattice calculation 
[T7] and chiral perturbation theory (TUl HHI (up to unresolved order ambiguities 
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Pni) and we retain it in the following. Table Ogives the final result for \Vus\, us- 
ing the KTeV slopes and the KLOE lifetime for Kl- For \Vud\ we use the value 
0.9740±0.0005 of Czarnecki, Marciano and Sirlin jTHj. The first error for each entry 



Table 1 : | Vus \ from kaon semileptonic decay. 



Source 


\Vus\ 


KTeV, Kl. e3,/i3 
KLOE, Kl. e3,/i3 
KLOE, Ks. e3 
E865, K+. e3 


0.2258±0.0013±0.0020 
0.2248±0.0015±0.0020 
0.2254±0.0017±0.0020 
0.2288±0.0021±0.0020 


Average 


0.2259±0.0022 


\Vud\ and unitarity 


0.2265±0.0022 



is due to experimental uncertainties on branching ratios, lifetime and form factor 
parameters. The common error of 0.0020 is due to the error on /+(0). Charged 
kaons and neutral kaons corrections are from ^ QUI ^j. The average is obtained 
by weighing according the first error. The x^-value is 2.5, for 3 dof, mostly due 
to the K~^ result which is 1.5a above the average. The error is the quadrature of 
0.0020, from /+(0), and 0.0008 from experiment. Unitarity is quite well satisfied: 
A = 1 - iKdP - |KsP=0.0003±0.0014, consistent with zero to 0.2a. The contribu- 
tions to the error on A from the two terms are equal. Fig. 4 gives a different view 
of the K'^ results. 
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Figure 3: The product /. 



Vus I for K^ decays compared with unitarity expectation. 



This way of presenting the data has the advantage of showing more clearly their 
consistency, without discrepancies being hidden by the large error on /:f'^''(0). The 
error /:f'^''(0) is now reflected in the width of the "unitarity" band. 

We do not include the result presented by Kleinknecht [201 at HQ&L2004 
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on June 2, because the Kf,3 branching ratio is obtained using the old NA31 value 
for BR{Kl 7r°7r°7r°) [21] • More recently [22], the NA48 presentation in Beijing 
at ICHEP04, using a judicious mixture of NA31 and KTeV values for the Stt*^ BR, 
gives, from the same data, the new preliminary value BR(i^e3)=0.4010±0.0045. 

8 \Vus\ from K fiu 

Marciano [22] has recently pointed out how the new precise lattice [21] results for 
the pseudoscalar meson decay constants, allow the determination of \Vus\ from the 
purely leptonic decays to muon-neutrino. In particular it is possible to relate the 
ratios |Kis|/|Kid| and T{K ^ iiv)/V{ti fiu). In the notation of eq. |21he gives: 



r{K ^ /x^(7)) 

r(7r /iz/(7)) 



ud 



f 

X ^ X /(masses) x {I + 6e 



P 

J TT 



V P fQl 
- 0.05271 ±0.0151 ^ ' 



\Vus\ = 0.2236 ± 0.0036 

In eqs. ^ fx/ fn = 1.201 ± 0.008 ± 0.015 comes from [23- The decay rates are from 
[7j and |V^rf|=0.9740±0.0005. From the previous discussion it is clear that the result 
for \ Vus\ is in good agreement with unitarity and also with the results quoted before. 
It is interesting to note that this last \Vus\ value is low while the one obtained from 
Ti^e^v is high. It is quite possible that the same cause is responsible for both 
effects. If the PDG K^2 "fit" value is low, because of old inaccurate measurements 
and/or improper accounting for radiation, than the value of \Vus\ will of course be 
low. As already pointed out in sec. HJ E865 has to use the PDG fit for the 
branching fraction in the 7r2, /i3 and 713 modes. If BR(iir^2) is low, then because of 
the fit 7r2+/i3+7r3 is correspondingly high, making the E865 final result also high 
by approximately Vl-5 as much, when expressed in per cent. KLOE has at present 
some preliminary indications that this might indeed be the case. 

9 Conclusion 

It is very gratifying to see the convergence of the new experimental results on kaon 
semileptonic decays. The CKM matrix appears now to be unitary and we can also 
use T^K"^ — > /i^z/(z/)(7)) and reach the same conclusion. While the previously 
claimed violation of unitarity is certainly not there, it is however disappointing that 
\Vus\ is still only known to about 1%. Agreement between different experiments is 
also still not as good as one could wish, but work is continuing. It is quite conceivable 
that in one or two years uncertainties due to lifetime, form factor shape parameters 
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and the branching ratios will be brought to the 0.1% level. We can only hope that 
a similar improvement will come about for /+(0). r decays do not agree with 
the direct observation in kaon decays. This is similar to r vs e"'"e~^7r+7r~ |2bj . 
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